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Abstract 

We extend the recent approach of M. Jimbo, K. Miki, T. Miwa, and A. Nakayashiki 
to derive an integral formula for the N-point correlation functions of arbitrary local 
operators of the antiferromagnetic spin-1 XXZ model. For this, we realize the quantum 
affine symmetry algebra Uq{su{2)2) of level 2 and its corresponding type I vertex oper- 
ators in terms of a deformed bosonic field free of a background charge, and a deformed 
fermionic field. Up to GSO type projections, the Fock space is already irreducible and 
therefore no BRST projections are involved. This means that no screening charges 
with their Jackson integrals are required. Consequently, our N-point correlation func- 
tions are given in terms of usual classical integrals only, just as those derived by Jimbo 
et al in the case of the spin-1/2 XXZ model through the Frenkel-Jing bosonization of 
Ug{su{2)i). 
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1 Introduction 



Solving many-body interacting quantum (statistical-mechanical or field-theoretic) systems 
is still a challenging problem even in low space-time dimensions. Recently, in a remarkable 
paper, Davies et al. [|I| presented a new method for solving the antiferromagnetic spin- 1/2 
XXZ model directly in the thermodynamic limit. This contrasts with the more conventional 
Bethe Ansatz approach, in which the theory is formulated on a finite chain before taking 
the thermodjTiamic limit in order to make the Bethe equations more tractable. Davis et 
al's approach is appealing for two obvious reasons: firstly, because in addition to being able 
to diagonalise the transfer matrix, it is also possible to calculate the correlation functions, 
form factors, eigenstates and S-matrix of the model; and secondly because it highlights 
the role and makes extensive use of the quantum affine symmetry of the model which is 
present only in the thermodynamic limit. The role of this symmetry is analogous to that of 
infinite-dimensional symmetries in conformal field theory. 

The Hamiltonian of the spin 1/2 XXZ Heisenberg quantum spin chain is [|], ^ ^ ^ 

1 oo 

nxxz = -^i: + ^f<i + A^f^r+i), (1-1) 

i=—oo 

where A = (g -|- q^^)/2 is an isotropy parameter. Davis et al. consider this model in the 



A < — 1 antiferromagnetic regime. The Hamiltonian |TT^ acts formally on the infinite tensor 
product, 

■■■V0V0V0V--- (1.2) 

where is a two dimensional representation of Ug{su{2)i) (in what follows the subscript k in 
Uq{su{2)k) or U'g{su{2)k) stands for the level of these algebras). The latter is the symmetry 
algebra of |lTl| and is a subalgebra (by dropping the grading operator) of the full quantum 
affine algebra Uq{su{2)i). The central idea of the approach of Davies et al. is to replace this 
infinite tensor product by the f/g(su(2)o) module 

^A,M = Hom(V(/i),V(A)), (1.3) 

where V"(A) and A = Ao,Ai are the Uq{su(2)i) highest weight modules and highest weights 
respectively. Following Frenkel and Reshetikhin Davies et al. realised this homomor- 
phism by introducing the following vertex operators (VOs) that intertwine the Uq{su{2)i) 



modules 

^'/{z)■.V{\)^V{^^)®V{z). (1.4) 

Here z is a spectral parameter and V\z) is the 'evaluation representation' of f/q(sM(2)o), and 
is isomorphic ioV ® C[z^ z^^]- These VOs 1.4 are the type I VOs of ref. [|^. By introducing 
a variant of these type I VOs, i.e., 

^^^y{z):V{\)®V{z)^V{^,), (1.5) 

Jimbo et al. were able to re-formulate the action of local operators of the vertex model 



(such as the spin local operator ol) on the infinite product |1.2| in terms of their action 
on J^A,/x- In this manner, and by bosonizing the VOs, they were able to write down an 
integral formula for N-point correlation functions of local operators in the 6 vertex model. 
Evaluating this integral for N=l, enabled them to triumphantly reproduce Baxter's result 
for the staggered polarization of the 6- vertex model 

Generalisations of the 6-vertex model, corresponding to spin-S" quantum chains, were 
introduced by Zamolodchikov and Fatteev and Kulish and Reshitikhin [IC]. The the- 
ory associated with 5 = 1 is a 19- vertex model |^, jTOf. The Hamiltonian is that of the 
antiferromagnetic XXZ model P, |3|, p| 



n = j \s^■s,^,- {s, ■ s,^,y + -iq- q-'Yis! ■ s^, - {s^ ■ s^.f + 2{sm (i.e) 

— {q + q ^ — 2) [(5*^5^ ■ S'^!,.]^ + Sf ■ Sjj^i)Sl ■ Sf^^^ + SI ■ S^_^_i{Si ■ S^_^_i + Sf ■ Sf^i)]^ , 

where S^, and are (3 x 3) spin-1 matrices, and J > 0. This Hamiltonian is symmetric 
under [/^(sm(2)2) [T^. The theory is massless for q = exp{i9), < 9 < n [l^, |14|, and for 
g = 1 it reduces to the XXX model discussed for example in [^, |15|. We consider it instead 
in the massive region with q real and — 1 < g < 1. The general spin-S* Hamiltonian is a 
complicated polynomial in the 25" + 1 dimensional spin matrices, and is given in the paper 
by Sogo ig. 

Idzumi et al. |]T2[ have extended the work of Davis et al. to general level k. The 
associated spin/vertex models are the generalised XXZ models with S = k/2 . The step 
of bosonizing the currents and vertex operators is however more complicated for k > 1. 



Firstly, one conventionally uses a deformation of the Wakimoto bosonization [116], O, IT8|, 



20] , pi] , - which requires three bosonic fields unhke the Frenkel-Jing bosonization for the 
case k = 1 [^], which requires only one. Secondly, and more seriously when calculating 
correlation functions it is necessary to take a trace over a Uq{su{2)k) irreducible highest 
weight module. For A; > 1, the bosonic Fock space can no longer be equated with any 
Uq{su{2)k) irreducible highest weight module, and one must therefore take its non-trivial 
BRST cohomology structure into account. This BRST structure can be studied by means of 



screening charges |2^, but the latter are given in terms of Jackson integrals, and a formula 
for the N-point correlation function derived using this approach would be highly complicated. 
In this paper, we realise the Uq{su{2)2) algebra and its type I vertex operators in terms of 
one boson and one fermion. The major advantage of this realisation is that now, up to a 
GSO-like projection (though not exactly the familiar one), the Fock space is irreducible as 
for k = 1. Exploiting this simplicity, we are able to produce an integral formula for the 
N-point correlation function of local operators of the spin-1 XXZ chain. 

The paper is organized as follows. In section 2, we fix the notation and make this paper 
self contained by reviewing the basic properties of the Uq{su{2)k) quantum affine algebra. For 
the purpose of free field realization we write this algebra in an operator product expansion 
form that we refer to as the Uq{su{2)k) quantum current algebra. In section 3, we use this 
quantum current algebra to realize Uq{su{2)2) and its associated type I vertex operators in 
terms of one boson and one fermion. Then we define the Fock space built from the modes 
of these two fields, and describe how, with the aid of a GSO-like projector, it is related to 
the highest weight representations of Uq{su{2)2). In section 4, we describe how the local 
operators of the spin-1 XXZ model act in JF, the space of excitation over the ground state. 
Here we introduce the variants of type I vertex operators. In section 5, we use the realization 
of section 3 to write an integral formula for the N-point correlation functions of arbitrary 
local operators. All integrals involved here are the usual classical ones and not Jackson 
integrals. Finally, section 6 is devoted to our conclusions. 
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2 The Uq{su{2)k) Quantum Current Algebra 

In this section we review some basic properties of the Uq{su{2)k) quantum affine algebra at 
level fc, which will make the subsequent sections more transparent. This algebra reads in the 
Chevalley basis as |2^, |2^, P7| 



Um'=q-^U, Uf^tr^ = q^f^, i^j, (2.7) 

[Ci, fj\ = SiJ q_q-l ; 

q'^aq-'^ = q^'-°ei, q'^fiq^'^ = q~^''°fi q\q~'^ = U, 

where {cj, fi, i,j = 0, 1} are the usual Chevalley generators and d is the grading operator. 
Its main feature is that it is an associative Hopf algebra with comultiplication 

A(ej) = Cj (g) 1 + tj (g) Ci, 

A(/.) = fi^ti' + l^fi, 

A{ti) = ti^ti, 

A{q'^) = q'^^q'^, 

and antipode 

a{ei) = -t-^e„ a(/i) = -/A, a(t,) = t;\ a{q'') = q-\ ^ = 0, 1. (2.9) 

When d is omitted, the above algebra is referred to as U'q{su{2)k)- For many practical 
purposes, it is convenient to use the Drinfeld realization of |2.7| . This is constructed from 
the following redefinitions of the Chevalley generators {e,, fi, tj}: 



to 




ti 


= 


eo 




ei 


= Eq, 


fo 




/i 


= Eo. 



(2.10) 
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Using |2.1CI| and |2.7| , one recursively generates the following algebra, which is known as the 
Drinfeld realization of Uq{su(2)k)'- 



[Hn, H„ 



[2n] ^nk_^-nk - 



2n q—q 

[Hn, Et] = ±V2^-^^^^Et^^, n ^ 0, 



[2.11] 



±V2Ho „±d ^±1/2 



This is an algebra generated by the Drinfeld operators {E^ (n G Z), Hm ijn G Z^o)? q 
where 7=*=^/^ commute with all the generators, and as usual [x] = (g^ — q^^)/{q — q^^)- ipn 
and ipn are the modes of fields ipi^) and ip{z) defined by 

E V'n^-" = q^""" exp{v^(g - q~^) E H^z-^}, 



,q 



n>0 



n>0 



E y^n;^-" = g-^""" exp{-v^(g - q-^) E i^n^^""}- 



2.12 



n<0 



n<0 



Using the basic identifications p.lO] one can re-express the comultiphcation |2]^ in terms of 

as 



the Drinfeld generators 2.11 



A{E+) = E+0 7^=" + ^^knqV2Ho 0E+ + Y.7=o 7^("+^'^/V„-i ® ^Hn^i)E^ mod ® A^^ 
ME^J = El 



'km 



7 + + Ello' I'-^V-m+i ® 7^^'"™^^i mod Ar„ ® Ar_2 



)7 



7'='" ® + ® g^^" + Eltl^ i^^'^-^^E^ ® 7-'=(™-^)/2^™_i mod iV! ® 



A(ff„) =Hrr,® 7'=™/2 ^ ^3fcm/2 ^ Hiod A^^ ® A^+, 

A{H^rn) = H_rn ® ^i'^^'"'!'^ + -^-^'^1'^ ® mod ® K 
A(7^2) = 7=*=" 



(2.13) 

where m > 0, n > 0, and N± and are left Q(g)[7^, ■i/'m, y9_„; m,n E Z>o] modules 
generated by {-E^j ^ £ Z} and m, n G Z} respectively 0, |2^. The main virtue of 

this comultiphcation is that it will make the derivation of the intertwining properties of the 



vertex operators explicit. These intertwining relations in turn allow the free field realizations 
(in terms of either free bosons, fermions, or ghosts) of the vertex operators. However for 
the purpose of the free field realization, it is convenient to rewrite the algebra in terms 
of operator product expansions (OPEs), that is, as a quantum current algebra (QCA). The 
Ug{su{2)k) QCA then reads [l^, HI, HI H 



ilj{z).(f{w) 

!f{z).E^{w) 
E+{z).E-{w) 
E^(z).E^(w) 



\'-'"iZ\i'-'"Cl'JM w).ij{z) 

±2(^-i»'?'F(2+fc/2)^ 



^-E^{w).iIj{z) 



(2.14) 



w{q 
{zq±^ 



z — wq 



z—wq" z—wq 

^E^(w).E^(z). 



I , \z\ > \wq 



±k\ 



Here the quantum currents ip{z) and y^{z) are given by 2.12 , whereas E^{z) are the following 
generating functions in terms of the Drinfeld generators: 



E^(z) 



+00 



(2.15) 



The free field realization of the Uq{su{2)k) QCA therefore translates into solving the 
relations |2.14| in terms of free fields. For example, the bosonization (i.e., the free field 
realization in terms of bosons) of Uq{su{2)k) has recently been studied intensively. For k = 1 
this is known as the Frenkel-Jing bosonization and requires only a single deformed boson 
field [pi]. For general k, there are many bosonizations available in the literature. We shall 
refer to these as q- deformations of the Wakimoto bosonization |l^, |l^, |18|, [19|, |2^, ^ |3^ . 
(See [0 for a detailed discussion of the equivalence of these different bosonizations.) Very 



recently, the bosonization of Uq{su{n)k) has been constructed in Ref. [PTI], and the difference 



realizations of Uq{su{n)) and general quantum Lie algebras have been achieved in Refs. 
and [B3| respectively. 



2.1 Intertwining relations of the vertex operators 

The vertex operators relevant to this discussion are the type I intertwiners of Refs [|I], [1^ . 
They are defined as maps between Uq{su{2)k) modules in the following way: 

$^'^'(2) : Vi\) ^ Vifx)®V^iz). (2.16) 
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Here ^(A) are Ug{su{2)k) left highest weight modules, with {A = Aj = (A; — i)AQ + iAi, i = 
0, . . . ,k} and {Aq, Ai} denoting the sets of Uq{su{2)k) dominant highest weights and funda- 
mental weights respectively. V^{z) {0 < j < k/2) is the spin j 'evaluation representation' of 
?7g(sM(2)o). It is isomorphic to ® C[z, z^^], where is the ?7g(s-u(2)o) 2j + 1 dimensional 
representation, with the basis {f^, ~ j < m < j} such that: 

where it is understood that vf^ = if \m\ > j. V^{z) is equipped with the following 



U'{su{2)o) module structure ||T2 





§2" 


= [j + € 








= [j - m]via+l € 








= [j - "^]^m+l S 








= [j + ^Vm-l ^ 








= q-^"'vi®z'', 






§2" 


= q^"'vi^®z''. 





(2.18) 



In terms of the Drinfeld realization this becomes 

7=^^/^t;4 ® = vi^® z^, 

E+vi®Z^ = g2n(l-m)[_^- ^/+n^ (2.19) 

Hnvl^ ®z^ = ^{[2nj] - + g-") [n{j + m)]}t;^ (g) z^+". 

Let V^* be the evaluation representation dual to and endowed with the left U'^^{su{2)q) 
module structure through the action of the antipode |2.9| , which is an anti-automorphism of 
this algebra, i.e., 

< xu,v >=< u,a{x)v >, X e Uq{su{2)o), ueV^*, v e . (2.20) 

Then it can easily be shown that the following two evaluation representations are isomorphic 
to each other |jl2|: 

C : ^4-. ^ Vr, (2.21) 



where 



Qj _ (^_'^y+mq-{j+m){j-m-l) 



2i 

j'+m 



, -j <m< j. 



(2.22) 



Here {vj^, —j<m<j} is the basis of V^*, which is dual to , and the notation 
defines the q-analogue of the binomial coefficient as 



(2.23) 



By definition the vertex operators (z) obey the intertwining condition [|], |I], 



(2.24) 



Here x denotes a generator in the Drinfeld realization, which is appropriate for explicitly 
constructing the vertex operators in terms of free fields. It is also convenient to define 
components of these vertex operators through 



m=-j 



(2.25) 



where the normalisation function {z) is given by 



_^^fc+2^|(A(A2,)+A(A)-A{M))^ 



with 



A(A 



2i. 



J(J + 1) 



(2.26) 



(2.27) 



k + 2 

{z) are referred to as the bare vertex operators. Using the above relation |2.24| , the 
arrive at the following commutation relations: 



comultiphcation |2.13| , and the fact that N^vLj = N_v-j = 0, N±vl^ G F[z,z ^]vi^^i, we 



[E+{w),(l>]iz)] 



0, 



[i/.,0^:(z)] = Jv/2{5„,o + gW/c+2)+|n|fc/2)l2M^n|^i(^)^ 



(2.28) 



'^^AZ), Eq ]g2j . . . , Eq ]q2(m + l) . 



In [^.28| there are (j — m) quantum commutators, where the quantum commutator [A, Sjqo; 
is defined by 

[A,B]q. = AB-q''BA. (2.29) 
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3 Realization of Uq{su{2)2) in Terms of One Boson and 
One Fermion 

When k = 2, one needs one boson, deformed in two different ways asx^(-z), and one deformed 
fermion "^{z) in order to realize the quantum current algebra |2.14| . The quantum currents 
are given in terms of these fields as 

= g2ao exp ( 2(g - q-') E anZ-A , 

\ n>0 / 

(p{z) = exp{ix^{zq~^) - ix^{zq)) (3.30) 
= exp ( -2(g - q-') E a^z-A , 

\ n<0 ) 

E^{z) = 72^(2) exp (±ix^(z)) , 



where 



[2n] ^„ \2n] 

n>0 ^ ' n<0 ^ ' 



^(z) = E Kz-^-^/\ 
rez+1/2 



(3.31) 



The bosonic modes {a, a„} of x^{^) and the fermionic modes br in the Neuveu Schwartz sector 
of ^E'(-z) satisfy the following Heisenberg algebra and anticommutation relations respectively: 



[a,ao] = i, (3.32) 

All the other commutation and anticommutation relations are trivial. x^{^) and "^{z) have 
the following simple OPEs: 

X^iz).x'^{w) = -\n{z - w)+ : x^iz).x'^{w) > 

X^(z).X^(w) = —\n{z — wq^'^)+:x'^{z)-X'^{'^)-: kl > (3.33) 
^(z).^(w) = ^7 irr^ ^+ ■.^(z)^(w) \z\>\q^^w\. 

It is clear from the above expressions that x^(^) (^(^)) are two different deformations of a 
usual real bosonic field (real fermionic field) and reduce to it in the limit q —>■ 1. The normal 
ordering symbol : ■ ■ ■ : introduced in |3.33| means that the annihilation modes {a„>o, &r>i/2} 
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are moved to the right of the creation modes {a„<o, &r<-i/2} and the shift mode a. In these, 
and in all other expressions in this paper, operators and products of operators defined at 
the same point z are understood to be normal ordered. The relation |3.30| is the realization 



in terms of one boson and one fermion of Uq{su{2)2) that we will be using from now on. 

When k = 2 the vertex operators can also be realized in terms of the same fermionic 
field "^{z) plus one deformed bosonic field ^(z) given in terms of the same modes {an, a} as 
To be more specific, let ^{z) be 

i{z) = a-^ao In {-zq') + 2z ^ ^ ^9 1 +^'^ To 1 ' (^.34) 

with 

= -Hq\wq^ - z))+ : az).^{w) :, \z\ > \q^wl (3.35) 

then the bare vertex operators are realized as follows: 

(f)i{z) = exp{i^{z)), 
Mz) = [Mz),E^]g. = V2§ti- Mz)E-{uj) : Io{uj,z)^{uj), 
= [M^),E,]/[2] = 2f ^ : M^)E-{u^)E-{v) : x 
Io{uj,z){n^{r],uj,z)'^{uj).'^{r]) + rl{v,oJ,z)^{7]).'^{uj)). 



(3.36) 



Here 



^o(^,^) = .,2(i_/,/4"i_V2^/,) , \zq<'\ < \w\ < \zq^\ 



The denominators are left in the form (1 — a) in order that the domain of convergence of 
the OPE sums \a\ < 1 is apparent. The w and rj integration contours must be chosen to 
lie within the domains indicated. This point will be discussed more in Section 5. Here the 
basic OPE's needed to derive ( |3.3(j| ) are given by 



h(z).E~(uj) =: M^)E-(ij): 



-1 



{zq'^){l~q-'^uj/z) ' 



E-iuj).Mz) =: E-{u)U^) : ^^^^^ (3.38) 
E''{uoi).E~{uo2) =: E^{uJi)E^{uo2) : {uJi - g^co'2), 
where E~{uS) indicates the purely bosonic part of E~[iS). In order to check our vertex 
operators ^.3(j| , we have computed two matrix elements involving the vertex operators 

•^Jr^''' (^) = (-2g')'/'(0+i(^) ® vl, + M^) ® + ® ^-1), t = 0,2. (3.39) 

10 



After performing the various double integrals, we find 

(3.40) 



The above results coincide exactly with those given in formula 3.27 of Ref. [|12| which were 
obtained by solving the quantum Knizhnik-Zamolodchikov equation. 

3.1 Fock spaces 



From the realization |3.3(J| and ^.36| it is clear that the representation of the Uq{su{2)2) are 
the Fock modules 

F{n) = F_^\n >, 
FJn > = 0, 



(3.41) 



where F^: are free Q(g) modules generated by {a=pm, ^t*-) ^ Z^o,r G Z>o + 1/2}. The 
states \n >, labelled by the integers n, are defined by 

|n >= exp(ma)|0 >, (3.42) 

where |0 > is the 'in' vacuum and it is annihilated by {am,m > 0;6,.,r > 1/2}. The dual 
Fock spaces are defined through at„ = a_„, af = a, and b]^ = b^r- F{n) is not irreducible 
but splits into two irreducible Fock modules F^{n) and F'^{n), which are isomorphic to the 
representations V{Xo) and V{X2) of Uq{su{2)2) respectively. The latter can be obtained from 
F{n) through the following GSO (though not exactly the usual GSO) projectors: 

P± = (1 ± exp {-27rid))/2 (3.43) 



as 



F°(n) = P+F(n) 
F'^(n) = P^F(n) 



(3.44) 



and where 



77,^ T \'2iT'\ 

- ^ = 4 E uT^a-nan + 2 E ITT^-r^- + («o) 72- (3.45) 

n>0 ^^'^J r>l/2 ^^'^J 

Moreover, F^{n) and F'^{n) are highest weight representations of Uq{su{2)2), with highest 
weights |0 > and |1 > respectively. The construction of the Fock module F^{n), which 
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is isomorphic to the representation V{Xi), is shghtly more comphcated but as it is 
unnecessary for our purposes we will not discuss it here. 



4 Correlation Functions of Local Operators 

A 'local operator' Ljv is one that acts naturally on the N tensor product V ^ V ^ . . . ^ V 
{V = V^l'^') corresponding to adjacent sites of a quantum spin chain. In order to interpret 
its action instead in JF = End{y[\i)\ which is the space of excitations above the ground 
state labelled by Aj, it is necessary to introduce another set of type I vertex operators 
denoted by ^^(Aj)y(-2) HI) Hi [H]- These vertex operators intertwine Uq{su{2)k) modules in 
the following way: 

^^(A.)y(^) ■■ ViaiX,)) ® Viz) ^ \/(A,). (4.46) 
Here \i = [k — i)Ao + iAi and cr(Ai) = iAo + [k — i)Ai. This VO has k + 1 components 
^aix.),v,J^) defined by 

^a{x,),v,J^)\^ >= ^a{x,),vi^)i\u > ®^m) = {idvix^® < v^, ■ >)^^^(^*^ iz)\u >, (4.47) 
where \u >G V{a{\i)) and 

^'(a!;(^) = <A.)(^^nA.) ® C)^":Uzq-'). (4.48) 

This implies that 

^t'ix.ivj^) = <AoC!e$i(^),_^(^g-'), (4.49) 
where ^^^^^^{z) = J2 ^t^{Xi)mi^) ® Vm and o.^\xi) ^ normalization constant, which we 



fix by ^ 



^aiX,),V,k/2-M^i^^) >= ^^^(X,)A^)iWi^^) > ^^k/2-^) = |A, > + ... (4.50) 

Using in addition the normalization 

^'(kV^/2(^)k(A0>=|A.>+... (4.51) 

we arrive then at 

1 



r 

^i-k/2 
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This in turn leads to 



C 



k/2 



~^l\x.)y,M = ^^Hi),-mi^^~')^ -k/2 <m< k/2. (4.53) 
If we set /c = 2 then cr(Ai) = X2-i and we get exphcitly 

^i^A^) = f.mi-^'izq-^ (4.54) 

where /«(+!) = /j(0) = — g*"^/[2] and /«(— 1) = g*. Here z = 0, 2 because only the 

ground states labelled by Aq and A2 |]12| will be considered in this paper. These ground 
states correspond to the two antiferromagnetic spin configurations, (■■■ + 1,— 1,+1,— I---) 
and (■ ■ ■ — 1, +1, —1, +1 ■ ■ ■) respectively. These are the ground states for which the staggered 
polarisation for the spin-1 XXZ model would be non-zero [|], 0. 



From the above definitions one can show that 

^x,.M^) ° ^xr'^i^) = 9xMv(x.). (4.55) 

where are scalar functions which can be obtained from the two-point matrix elements 
|3.40 and relation 4.54 . We find then 

gx, = —^. z = 0,2. (4.56) 
1 — 



The action of L^v on JF can be specified by using the VOs ^^/{z) and <l^y(2) 0, HI 
the following way: 



= i9x.r' (4.57) 

(-i)^5::iv^(-2) • ■ ■ ^i:::y^ i^N){^dy,^,,^ ® (-n) ■ ■ ■ $5-^ (.,)$5r -^(.o 

Here z = or 2 labels the choice of ground state (or equivalently the Uq{su{2)2) module over 
which the trace is to be taken), zi, - ■ ■ ,zn are local spectral parameters, and the subscripts 
on the As are understood as modulo 4. Then the correlation function of this operator is 
given by: 

<L>= iJe^::ip:^:::i{zi, ^2, ■ ■ ■ , zn\i), (4.58) 

where 



P:^:::,'{z,, Z2,---, zn\i) = (1 - qYx (4.59) 

\Xi ( \ ^Ai+2{iV-l) / \^\+2Ny , 
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and p = 2d + 0, |T2| . This trace may be rewritten purely in terms of the bare vertex 

as 



operators of |3.36 



N 



1=1 

T"<--'^'^ ''{ziq-\ ZNq~\ zn,--, Zi\x, y) + 



2ao' 



(Tr(a;-V"") + (-l)'/'Tr(x-'^?/ 
where x = g^, X"*^/^ = — x^^^, ?/ = g^^ and 



The traces in |4.6CI| and [4.61| are taken over the complete Fock space of Section pT . 



(4.60) 



(4.6i; 



5 Evaluation of the Trace 



The trace 4.61 splits naturally into a product of contributions from the bosonic non-zero 
modes, fermionic modes and bosonic zero modes |3^. The evaluation of these three terms 
will be dealt with separately. 



5.1 The trace over bosonic non-zero modes 



From the expressions |3.36| for (j)i{z) , (f)Q{z) and (j)^i{z), it is apparent that the complete 
contribution of bosonic non-zero modes to [4.61 



IS 



where = 4 i^a-nCin, and 

n>0 ' ' 



Fe. = l 
Fe,=0 

F.,=-i 



(5.62) 



i=l 



Mzi), 

: E~{uJi)(f)i{zi) :, 



(5.63) 
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The hat indicates the bosonic non-zero mode part, and lui and rji are integration variables. 
To evaluate ^.621 we use the trace reduction technique of Clavelli and Shapiro . The trace 
of an operator O which is a function of the bosonic non-zero modes a„, is given by 

Trfx-^^O) = 5.2M^^ (5.64) 
[x; X) 

oo 

where [s; x) = Y\ {1 — sx"), and O is the same operator expressed in terms of oscillators a„ 

n=0 

defined by 

an = z — '^ + c^n{n>0), an = an + ——^ {n < 0). (5.65) 
1 — x"- — 1 

Here an extra set of oscillators c„ have been introduced which have the same commutation 
relations as the a„, and commute with them. The technique for calculating < 0|(9|0 > is to 
completely normal order the operator O with respect to a„ and c„, and then to use Wick's 
theorem (for exponentials) in order to express it in terms of contractions between pairs of 
operators. When carrying out the integrations associated with (f)o{z) and 0_i(z), the contour 
must always be chosen to be consistent with the domain of convergence of the series that 
occur in the normal ordering calculation. Specifically this means that whenever 1/(1 — au) 
or {aw; x) appears on the right-hand side of an OPE the contour for integration over u must 
be chosen such that \au!\ < 1. 

After normal ordering with respect to a„ and c„, (f)i{z) and E~{uj) become 



n<0 



[2n] J \ ^0 il-x-)[2n] 



N /2 N / (fln^ "-Cna;")g "\ / ^ (a„^ " - c„(a;x)")g"\ 
E (tu) = (q x; x) exp 2 > — exp 2 > ; — . (5. 

Re-normal ordering with respect to a„ and Cn (indicated by o ■ ■ ■ o), gives 

(g^x; xY 



E-{uj)<Pi{z) : 



g{uj/z) 



: E {ui)E {u2)Mz) ■■ = o-.-o— — — - — — , 5.67 

g{uji z)g{uj2 z) 



where 

g{uj/z) = {q-^{uj/z)x;x){q^{z/uj)x;x), 
h{ui/uj2) = {q^{ui/uJ2)x]x){q^{uj2/^^i)x]x). 



(5.68) 
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Normal ordering in pairs gives 



where 



0(j)l{Zl) O 0(f)i{z2)o = 0---oGl{zi,Z2), 
0(j)i(^z) O oE^ {lj)o = O ■ ■ ■ O G2{Z,U), 

o£'^(C(j) o o0]^(z)o = o ■ ■ ■ o Gsi^uj, z), 
oE^ {uoi) o oE~ {002)0 = o ■■■ o 04(001,002), 

Gi{zi,Z2) = {q'^{z2/zi);x){q'^{zi/z2)x;x), 

G3{UJ,Z) 



(5.69) 



1 



(5.70) 



{q^ {z /Lu);x){q^'^ {lj / z)x;x) ' 

G4{uJi,U2) = {q^{u2/uJi);x){q^{uJi/uJ2)x;x). 
To label the fields we use the indices 

a, a' & A = {/|ei = 1, 0, — 1}, dim A = n, 

b, b'eB = {l\ei = 0,-l}, dim B = nB, (5-71) 

c, c' e C = {l\ei = -1}, dim C = nc- 

Then in terms of these functions the non-zero bosonic mode contribution 5.62 to the trace 
is 



4.61 



i=i \x;x) b gi^b/ Zb) c givj Zc) 

n GiiZa,Za>)Y[G2iZa,iOb)Y[G2iZa,Vc) (5.72) 

a<a' a<b a<c 

Yl Gs^Ub, Za) n Giiujb, ujb') n G^iuJb, r]c) 

b<a b<b' b<c 

Yl G^{ric, Za) n Giij]^, tOb) n G4iVc, Vc), 

c<a c<b c<c' 



where H = 11 etc. 

b b&B 



5.2 The fermionic trace 

Introducing the set of integers {/c} with 1^ = 1 or —1 for each c G C, which shall label both 
the order of the fermions and the associated l'^i(ric,iOc, z^) functions of |3.36| , the fermionic 



contribution to [4.61| becomes 

Trix-^'fUHb), (5.73) 



b 
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where —dj = 2 i^^-r^r, and 

r>l/2 ' 



Hb={ 



^{ub)'^irit), beC, h = l, 
^(r/fe)$(cjb), beC, lb = -l, 



Using Clavelli and Shapiro's formula for fermionic trace reduction we get 



Tr{x-'^fO) =< 0|O|0 > l[{l + x'). (5.74) 



r>-k 



Now O imphes br br, where, 



+d-r ir>h , br = br + -^-^ ir<-h. (5.75) 



Here dr are anticommuting copies of the br- Normal ordering with respect to br and dr gives 

^{Ujl).-^{UJ2) = 0---0 + A{UJI,UJ2), (5.76) 

where, 

A(.,.,) = i^2:(,- + ,--)i(^^lM^J^2A^, (5,77) 

^ r->i [i- + X ) 

' — 2 

This expression can be split into two parts, Ai(co'i, 0^2) and A_i(c<Ji, c<j2), where, 

^ r> 1 (^i + X ) 

' - 2 

Ai(c<Ji, c<j2), and A_i(to'i, ^2) are convergent in the strips |x| < \c^ijJ\IijJ-2\ < 1 and 
|x| < \q^'^uj2/uji\ < 1 respectively. For the purpose of integration Ai{ui,U2) may be re- 
expressed in terms of the usual 9 functions through the identity of Goddard and Waltz 



Aiicu„U2) = '^""^^^^ ' e2iO\T)eMr)03iui\r)/9MT), (5.79) 

where ui = ln(g^'c(j2/co'i)/(27ri) and r = ln{x) / (2TTi) . Writing the theta functions in terms of 
infinite products gives 

_ {x'^;x^f{x^;xf {-x^q^^uj2/uJi)]x){-x^q'^^uj2/^^i)~^;x) 
1[UJ,,UJ2) 2{q-^co^ - q^co2) (x(g2^a;2M); x) (x(g2^a;2M)-'; x) ' ^^'''''^ 
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Then |5.73| is given by 

Tr{x~'^fY[Hb) = < OlH^fclO > (5.81) 

b b 

where the vacuum-vacuum correlation function is given by Wick's theorem as a sum over all 
possible products of pair wise contractions of the fermionic fields. Each of these contractions 
is given in terms the fermionic "propagator" 

< ^(a;i)^(tU2) >= Ai{uJi, UJ2) + A_i{ui, UJ2). (5.82) 



5.3 The trace over zero modes 



The contribution of the bosonic zero modes to t4.61| is the factor 



i=l 



where —do = al/2, and in analogy to equations p.63| we have 

= :Eo-(u;.)0?(^.):, 
with the zero sub /superscript indicating zero modes. We define the functions 



(5.83) 



(5.84) 



0?(^l)0^(^2) 

Eo{uj)4>\z) 



:0?(^)Eo"(^):G'O(z,a;), 
:Eo'(^)0O(z):GO(u;,z), 



4 



G\{zi, Z2) 
Gl{z,uj) = {-zq')-\ 



Gl{u,z) 



Eq {uJi)Eo {LJ2) = : Eq {ui)Eq {UJ2) : Gl{uJi,uj2), Gl{uJi,U2) 



(5.85) 



After completely normal ordering all operators in ^.83| , and provided that = 0, without 

i=l,n 

which the trace vanishes, the factor left over is 



Tr x'^oy'^o nexp(ao ln(-z,g^)) nexp(-ao In(cUfe)) nexp(-ao Hr],)) 



(5.86) 
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Thus p.83| is given by 



i=l 



mez 



tor] 



a<a' a<b a<c 

n Giicob^za) n G'i{iUb,tob') n ^0(^^5,77,) 

6<a b<b' b<c 

n GO(r/e,Za) n GlivciOb) n GO(r7„r?,0, 

c<a c<6 c<c' 



where z = Y\Za, uj = JJoOb and ?7 = 11 

a b c 



(5.87) 



5.4 Collecting terms 

Defining Gi = GiG^^ (/ = 1, ■ ■ ■ , 4), the complete expression for 4.61 is given by 



n / ^^o(^., ^'') E {n / < oi n ^^io > 

1 



b 

n 



n 



h gi^b/zb) c givc/zc) 

Yl Gi{Za,Za')Y[G2{Za,UJb)Y[G2{Za,r]c) 

a<a' a<b a<c 

Yl Gs^Ub, Za) n Gi{ujb, t^b') n GiiuJb, Vc] 

b<a b<b' b<c 

Y[ G-i{r]c, Za) Yl GiiVc, tUfo) n G^iVc, Vc) 



c<b 



c<c' 



(5.88) 



where the first sum is over all sets of {Ic} with each Ic = ±1. The contour for each integral 



in this sum is chosen in the way described in section |5.1| . Each contributing integral will 
correspond to one choice of {Ic}, and a particular fermionic contraction. Using the Jacobi 
triple product identity p8[] , 



(5.89) 



meZ 



the last factor in |5.88| within [■ ■ ■] can be written as 



(x2;x)(-x2( ^ ' ); x){-x^ C - )" ; a;) 



our] 



(5.90) 
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Similarly, 



J2 (x'^'^V"") 

meZ 



(x; x) 

(x2 ; x){—x2y'^; x){—x^y~'^; x) 



(5.91) 



N-point correlation functions of local operators are then given in terms of 4.58, 4.60 and 



6 Conclusions 

In this paper, we have derived an integral formula for the N-point correlation functions of 
arbitrary local operators of the antiferromagnetic spin-1 XXZ model in the thermodynamic 
limit. In doing this, we have derived a one boson one fermion realization of Uq{su{2)2). 
This realization is more convenient than the q-deformation of the Wakimoto realization of 
this algebra, which is given in terms of three deformed bosons. This is because the latter 
realization leads to infinitely many Fock spaces that are highly reducible. To project out 
the irreducible subspaces, which are isomorphic to the highest weight representations of 
Uq{su{2)2), one has to study the cohomology structure of the Fock spaces by means of the 
screening charges, which act as BRST operators among them. These screening charges are 
non-local and given as Jackson integrals in terms of the screening currents. This means 
that had we used the q-deformation of the Wakimoto construction, the N-point correlation 
functions would have involved, in addition to the usual integrals, both Jackson integrals and 
infinite summations arising from the BRST projections. These two would be shortcomings 
are avoided in our one boson one fermion realization of Uq{su{2)2) because the Fock spaces 
are already irreducible (up to simple GSO-like projections as explained in section 3). 

It is still a challenging problem to try to explicitly integrate, if not the N-point correlations 
functions, at least the 1-point correlation functions. Choosing the local operator as 
would lead, as for k = 1 0, to the staggered polarisation of the spin-1 XXZ model. Other 
quantities of physical interest, which might be extracted with this technique, are the local 
energy density, and form factors of the model. 
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